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In the nuclear reactor design, a code for automatically generated multi-temperature continuous-energy neu- 
tron cross section data library, which is called AMTND for short, was designed and developed to meet the need 
of the reactor core design coupled with thermal-hydraulic design. The code can provide a point-wise cross- 
section at any temperature for a Monte Carlo neutron transport program, such as MCNP. In ensuring that the 
nuclear data produced by AMTND meets the testing of critical benchmark experiments, the time-consumed by 
the nuclear data generating of AMTND compared with NJOY’s was carried out and the result shows the code’s 
excellence. In order to test the accuracy of the code, the Doppler coefficient test benchmark was also carried 


out and the results verified the code preliminarily. 
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I. INTRODUCTION 


Nuclear data libraries are the basic of neutronic design 
for reactors. In the core design, various reactor parameters, 
for example, reactivity coefficients of the reactor transition 
from a cold state under zero power to a hot state under zero 
power and to a hot state under full power, needed to be simu- 
lated. The core design will also be coupled and iterated mutu- 
ally with thermal-hydraulic design. In the coupling process, 
each set of thermal-hydraulic parameters should correspond 
to a set of reactor design parameters. A variety of nuclear 
databases for real-time temperature needs to be generated in 
calculating these parameters of the reactor core. 

In this paper, based on reactor design studies at home and 
abroad and the calculations in the actual operation process, 
an Automatically generated Multi-Temperature continuous- 
energy Neutron cross section Data library code, called 
AMTND for short, was designed and developed. 


Il. THE DESIGN OF AMTND 
A. Background 


Using the Monte Carlo neutron transport code (MCNP) is 
widely recognized as more precise and mature than any other 
neutron transport program code in reactor design and anal- 
ysis. The point-wise continuous-energy cross section data 
in ACE format is suitable for use by the MCNP family of 
Monte-Carlo codes and other codes that can read the same 
format library. 
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When using MCNP to calculate temperature-related pa- 
rameters in the design of a rector core, nuclear data from cor- 
responding temperatures needed to be produced to ensure its 
accuracy. That was quite a waste of time. A nuclear data 
processing package SIGACE [1] was released by IAEA. It 
could broaden the ACE format data at merely 300 K to higher 
temperatures based on FENDL2.1 nuclear data. But it would 
produce errors when dealing with some heavy nuclei, such as 
the fallacious values of the nuclear data’s location parameter 
JXS array [2]. 


B. Doppler-broadening 


In the neutron calculations, the effective cross section for 
a material at temperature T is defined as the cross section 
that gives the same reaction rate for stationary target nuclei as 
the real cross section gives for moving nuclei. In many cases 
of interest, the target motion is isotropic and the distribution 
of velocities can be described by the Maxwell-Boltzmann 
function [2]. Therefore, the standard form of the Doppler- 
broadened cross section is as follows: 


a I dVo(V)V? oe = ae 
V° Jo 
(1) 


where, a = M/(2kT), V = |v — v'|, k is Boltzmann’s con- 
stant, M is the target mass, v is the velocity of the incident 
particles, v’ is the velocity of the target, and o(T,v) is the 
cross section for nuclides. 


o(T,v) = 


C. The preparation of the MCNP input file 


Writing the information for radionuclide species and ra- 
dionuclide temperatures that needed in the MCNP neutron 
transport calculations in the MCNP input file without inter- 
fering with the normal operation of MCNP. For example: 
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M1 82208.83c 3.91128E-3 $ t=356.9 K 


This is the description of the material card of the MCNP input 
file. The original MCNP input file is on the left of the charac- 
ter ‘$’ and the comment is on the right. ‘82208’ stands for the 
radionuclide types of elements. The first two digits are for the 
number of protons and the last three digits are for the mass of 
the nuclides. ‘83c’ is for the database identifier. It can be 
defined according to the user’s own needs. In order to dis- 
tinguish a normal radionuclide from a transient radionuclide, 
the transient radionuclide was defaulted to ‘*1.c’ deliberately. 
‘t=356.9 K is the nuclides temperature. All of the letters are 
not case sensitive. 


D. The programming of AMTND 


This work completes the production of the database library 
and the neutron transport calculations in the process. The 
flow chart of AMTND is given in Fig. 1. 


1. The AMTND reads radionuclide species information, 
database identifiers, and radionuclide temperature from 
the material card of the MCNP input file. Meanwhile, 
AMTND will merge the identical information of all the 
nuclides automatically and store it in a temporary struc- 
ture array. 


2. AMTND sealed the radionuclide species information, 
radionuclide identification (ID) number, and radionu- 
clide evaluation database name (408 nuclides come 
from HENL-ADS) in a structure array. Using the ra- 
dionuclide species information obtained from MCNP 
input file and the radionuclide ID number, a radionu- 
clide evaluation database name could be found in the 
structure array above. 


3. Using the information on the radionuclide species, ra- 
dionuclide ID number, database identifiers, and the 
temperature obtained above, AMTND generated a 
Doppler-broaden input file automatically. 


4. AMTND calls for Doppler-broaden modules. It reads 
the Doppler input file and the evaluation database files 
and then generates ACE format files and corresponding 
index files. 


5. Merge the ACE format files and the corresponding in- 
dex files into the original database. The new database 
can be used by MCNP. 


Il. TESTING AND DISCUSSION 


The testing of the code was mainly reflected in the accu- 
racy of the data and the time consumed by generating the nu- 
clear data. There are three tests as follows. 
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Fig. 1. The flow chart of AMTND. 


A. Critical benchmark experiments 


Critical benchmarks testing mainly deals with some ther- 
mal and fast reactor data [3—8]. In this test, 10 fast reactor 
experiments, which were released by Organization for Eco- 
nomic Cooperation and Development (OECD), were selected 
from the “International Handbook of Evaluated Criticality 
Safety Benchmark Experiments”. The evaluations of the ra- 
dionuclide come from HENDL-ADS [9] and the temperature 
is 293.6 K. 

Based on criticality calculations using MCNP with the data 
produced by AMTND, the overall agreement of Keg between 
experimentally measured and computed shown in Table | was 
very good. Some of the deviation between the calculation 
results and the experiment values was less than that produced 
by NJOY. 


B. Time-consumed by generating nuclear data 


As it is shown in Table 2, AMTND has been greatly im- 
proved by NJOY99 in generating an actinide nuclei data li- 
brary. The average speedup ratio is 1.37 and the maximum 
speedup ratio is 1.61. 


C. SEFOR Doppler benchmark 


The SEFOR reactor [11] was designed to provide a 
Doppler measurement in an environment that is representa- 
tive of an operating LMFBR with respect to the neutron spec- 
trum, fuel temperature range, reactor composition, and the 
fuel microstructure. Standard fuel for SEFOR was mixed 
oxide (20% PuO,, 80% UO,), in which the Pu contained a 
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TABLE 1. The critical test of database produced by NJOY99 and AMTND 


Experiment ID NJOY-MCNP AMTND-MCNP Experimental values Deviation (%) 
U233-MET-FAST-001 0.99849 1.00019 1.0000 0.019% 
U233-MET-FAST-002-case1 0.99832 0.99921 1.0000 0.079% 
U233-MET-FAST-004-case1 1.00500 1.00417 1.0000 0.417% 
HEU-MET-FAST-001 1.00111 1.00028 1.0000 0.028% 
HEU-MET-FAST-027 1.00492 1.00394 1.0000 0.394% 
PU-MET-FAST-001 0.99981 0.99989 1.0000 0.011% 
PU-MET-FAST-005 1.00013 0.99780 1.0000 0.220% 
PU-MET-FAST-035 1.00532 1.00417 1.0000 0.417% 
SPEC-MET-FAST-008-case1 1.00381 1.00335 1.0026 0.335% 
SPEC-MET-FAST-001 1.00198 1.00216 1.0000 0.216% 


TABLE 2. The time consumed by the nuclide data produced by TABLE 3. SEFOR Doppler benchmark, Spherical Model Region 


NJOY99 and AMTND in Doppler broadening (s) composition (atom/barn-cm) 
NJOY99(t1)" AMTND(t2 J speedup ratio (tı /t2 " Material Region] Region2 Region3 Region4 
BBY 25.2 18.7 1.35 Fe 1.3574E-2 1.3886E-2 5.8932E-3 7.8587E-3 
235, 16.1 13.2 1.22 Cr 3.9574E-3  3.9511E-3 2.8913E-3 2.4623E-3 
2333y 7.9 7.6 1.04 Ni 2.0292E-3 2.3580E-3 3.0178E-2 1.3315E-3 
232Th 15.7 10.9 1.44 Na 1.6615E-2 6.8099E-3 5.4493E-3 1.3070E-3 
239p 17.5 11.3 1.55 Be — 3.6011E-3 1.8327E-5 - 
240Dy 12.9 8.0 1.61 O — 2.0991E-2 1.2597E—4 — 
— S Mo — 1.1999E-4 1.5605E-5 — 
aà ti is the time NJOY99 consumed; B-10 — 6.1100E-5 = 5.7684E-3 
P t9 is the time AMTND consumed; B-11 = 2.4600E—4 = 2.3100E—2 
€ The evaluation of radionuclide are ENDF/B-VII [10]; U_235 7 1.5374E-5 1.1724E-7 _ 
d CPU, Intel®Core™i5-24006. U_238 _ 6.9808E-3 53438E-5 _ 
Pu-239 — 1.5901E-3 — — 
Pu-240 — 1.4355E—4 — — 
Al — 1.4355E—4 — — 
C — — 2.2330E-3 6.5800E-3 


TABLE 4. The calculation results of SEFOR model 
kı k2 Tdk/dT 


NJOY 1.006 10+0.00055 1.000 77+0.000 56 —0.0076+0.0011 


a AMTND 1.007 090.000 57 1.001 810.000 58 —0.0075+0.0012 


0.0 10.688 49.75 79.185 92.666 


kz =neutron multiplication factor with the fuel at T>. 
Fig. 2. 1-D SEFOR Doppler benchmark spherical model (in cm). The radionuclide of the SEFOR model adopted is the 
evaluation documents of HENDL-ADS/MC [9]. The re- 
minimum amount of ?Pu (~8% of the Pu) and the U was sults are shown in Table 4. A correction factor for model- 
depleted in ?5U. A one-dimensional spherical model is ana- ing of —0.0005 was obtained due to the resonance hetero- 
lyzed in the present study. The 1-D SEFOR Doppler bench- geneity, subassembly heterogeneity and other factors [11]. 
mark spherical model is shown in Fig. 2 and the spherical The calculated value of the Doppler coefficients were 
model region composition is shown in Table 3. More details —0.0081+0.0011 using NJOY-MCNP and —0.0080+0.0012 
can be found in Ref. [11]. Using the Doppler Calculation using AMTND-MCNP. These results compared well with the 
Model, the isothermal Doppler coefficient is computed as de- experimental result of —0.00800.0010. 

scribed below: 


it b-i ae IV. CONCLUSION 


z = , (2) 
ie) EONO In this work, an automatically generated continuous- 
Tı = 677 K; average fuel temperature at zero power; energy neutron cross section data library code called 
Tə =1365 K; average fuel temperature at 20 MW; AMTND was designed and developed. The code was tested 
kı =neutron multiplication factor with the fuel at T3; by critical benchmark experiments, and was time-consuming 
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from generating the nuclear data and the isothermal Doppler 
coefficient test benchmark. The results preliminarily showed 
the accuracy and reliability of the code. More verification 
about the code will be needed in the future. 
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